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ABSTRACT 


A  liquid  crystal  thermographic  technique  has  been  devel- 
oped which  provides  an  excellent  means  of  obtaining  both 
qualitative  and  quantitative  heat  transfer  information  on 
heated  objects  placed  in  forced  convection  environments. 

Circumferential  variation  of  the  Nusselt  number  on  a  uni- 
formly heated  right  circular  cylinder  cooled  by  forced  con- 
vection was  obtained  for  Reynolds  numbers  varying  from  38,000 
to  148,000.   The  results  compare  within  the  experimental  un- 
certainty in  forward  stagnation  regions  with  the  theory  of 
Schuh.   Beyond  approximately  30°,  the  results  diverge  from 
theory  but  are  consistent  with  the  work  of  other  investiga- 
tors . 
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I.   INTRODUCTION 

The  purpose  of  this  investigation  was  to  develop  a  tech- 
nique for  the  visual  determination  of  both  qualitative  and 
quantitative  heat  transfer  and  fluid  flow  information  on 
heated  objects  placed  in  forced  convection  environments. 
Cholesteric  liquid  crystals,  a  material  that  exhibits  bril- 
liant changes  in  color  over  discrete,  reproducible  tempera- 
ture bands,  were  used  as  the  temperature  sensor  in  the 
technique.   The  liquid  crystals  allow  one  to  visually  observe 
select  isotherms  and,  further,  can  be  used  to  infer  the  loca- 
tion of  points  of  flow  separation  and  boundary  layer  reattach- 
ment.  The  liquid  crystals  also  give  a  dramatic  indication 
of  the  influence  of  turbulence  on  surface  temperature. 

The  technique  has  been  used  to  determine  the  circumferen- 
tial variation  of  the  Nusselt  number  on  a  uniformly  heated 
right  circular  cylinder  cooled  by  forced  convection  in  a 
cross  flow  of  air.   Reynolds  numbers  were  varied  from  38,000 
to  148,000  in  the  present  investigation.   This  allowed  the 
study  of  both  critical  and  subcritical  flows.   The  free 
stream  turbulence  intensity  was  approximately  0.5  -  0.7%. 

Data  were  obtained  on  a  4  inch  diameter  right  circular 
cylinder  constructed  from  .039  inch  thick  carbon  impregnated 
paper  that  had  a  resistivity  of  approximately  1  ohm-in.   The 
surface  of  the  cylinder  was  electrically  heated  by  passing  a 
known  current  longitudinally  through  the  paper.   The  inner 
hollow  space  formed  by  the  cylinder  was  firmly  packed  with 
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glass  wool  to  prevent  heat  losses  into  the  cylinder.   The 
glass  wool  also  reinforced  the  cylinder  and  aided  in  resisting 
deformation  due  to  the  outer  flow. 

The  results  obtained  in  the  present  investigation  com- 
pared within  the  estimated  experimental  uncertainty  in  the 
forward  stagnation  region  with  a  theory  proposed  by  Schuh  [1] . 
Beyond  approximately  30°,  the  experimental  values  rapidly 
diverge  from  the  theory.   This  is  consistent  with  results 
found  by  other  investigators  [2,3]  and  is  most  probably  ex- 
plained by  the  fact  that  the  pressure  distribution  used  in 
Schuh* s  theory  was  taken  to  be  the  ideal,  frictionless  dis- 
tribution. 

II.   BACKGROUND 

The  use  of  liquid  crystals  in  wind  tunnel  experiments  was 
first  investigated  by  Klein  [4]  in  1968.   Although  they  showed 
great  promise  for  obtaining  heat  transfer  rates,  several 
problems  existed.   The  degree  of  accuracy  was  limited  by  the 
influences  of  mechanical  shear,  ultraviolet  light,  and  chem- 
ical contamination. 

Some  time  after  the  completion  of  Klein's  study,  the 
National  Cash  Register  Company  (NCR)  developed  a  process  for 
encapsulating  liquid  crystals.   McElderry  [5]  used  these 
crystals  in  a  study  of  boundary  layer  transition  at  super- 
sonic speeds,  and  found  the  encapsulated  crystals  to  be 
relatively  insensitive  to  shear  and  contamination. 


Wirzburger  [6],  in  a  thesis  study,  used  liquid  crystals 
to  study  the  heat  transfer  characteristics  of  a  container 
stored  rocket  motor  placed  in  a  hostile  environment.   The 
temperature  distributions  produced  by  resistively  heated, 
cryogenic,  and  radio-frequency  surgical  probes  were  studied 
by  Groff,  Petrovic,  and  Katz  [7,8,9,10,11,  12]  using  liquid 
crystals. 

One  of  the  original  objectives  of  a  thesis  investigation 
by  Meyer  [3]  was  to  use  encapsulated  liquid  crystals  to  study 
the  surface  temperature  distribution  on  a  uniformly  heated 
cylinder  cooled  by  a  cross  flow  of  air.   This  objective  was 
never  accomplished.   The  temperature  distributions  in  his 
study  were  all  found  using  thermocouples.   The  present  inves- 
tigation is  a  continuation  of  the  work  initiated  by  Meyer. 

III.   EXPERIMENTAL  APPARATUS 

A.   LIQUID  CRYSTALS 

There  are  three  types  of  liquid  crystals;  smectic,  nem- 
atic,  and  cholesteric  [13]  .   Of  these  three,  only  the  choles- 
teric  type  was  used  in  the  present  investigation. 

Microencapsulated  liquid  crystals  provide  a  visual  dis- 
play of  temperature  when  applied  to  a  black  surface.   The 
liquid  crystal  surface  color  will  pass  through  the  visible 
color  spectrum  in  sequential  order  as  a  function  of  surface 
temperature.   The  temperature  at  which  a  given  color  appears 
is  a  function  of  the  cholesterol  ester  formulation.   By 
selecting  the  proper  formulation  the  entire  spectrum  (red  to 
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violet)  can  be  transversed  in  fraction  or  multiple  degrees 
with  an  accuracy  of  approximately  0.1°C  [14].   Fergason  [13, 
15-17],  in  a  series  of  articles,  presents  excellent  discus- 
sions on  the  chemistry,  varieties,  properties,  uses,  and 
limitations  of  liquid  crystals. 

The  process  of  microencapsulation  consistsof  enclosing 
the  liquid  crystals  in  20-30  micron  polyvinyl  alcohol  capsules. 
These  capsules  are  suspended  in  a  water  based  slurry.   This 
process  makes  the  liquid  crystals  relatively  insensitive  to 
mechanical  shear  and  chemical  contamination.   It  also  reduces 
the  sensitivity  to  viewing  angle  when  compared  to  raw  liquid 
crystals  [5] . 

Fourteen  different  liquid  crystal  temperature  formula- 
tions were  used  in  the  present  study.  A  temperature  range 
of  30°C  to  49°C  was  covered.  The  liquid  crystals  used  and 
their  calibration  points  are  listed  in  Table  I. 

The  liquid  crystals  were  calibrated  using  a  Rosemount 
constant  temperature  bath  following  the  general  procedure 
outlined  by  Petrovic  [8] . 

All  crystals  were  calibrated  on  a  piece  of  the  material 
to  which  they  would  be  applied  for  data  collection.   With 
the  exception  of  liquid  crystal  R-45  all  liquid  crystals  in 
Table  I  were  calibrated  on  a  piece  of  the  carbon  impregnated 
paper.   A  piece  of  tape  was  the  base  material  for  the  R-4  5 
calibration. 
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Table  I 
Liquid  Crystal  Calibration  Results 


Liquid  Crystal 
(NCR  Designation) 

S-30 

S-32 

R-33 


S-34 
S-36 
R-37 
S-38 
^-40 
R-41 
S-43 
S-45 
R-45 
R-49 


Color  Transition 

Red(°F) 

Green (° 

F) 

Blue(°F) 

85.9 

86.4 

87.5 

90.5 

91.6 

92.6 

92.8 

93.7 

95.6 

94.4 

95.2 

96.3 

97.5 

98.2 

.  99.4 

98.9 

99.7 

102.4 

101.4 

102.4 

103.4 

105.1 

106.1 

107.1 

106.1 

107.4 

109.8 

109.8 

111.0 

112.3 

113.2 

114.0 

115.1 

108.9 

109.9 

121.1 

121.7 

122.0 

124.2 
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B.   CYLINDER  DESIGNS  AND  CONSTRUCTION 

The  basic  experimental  procedure  was  first  developed 
using  the  cylinder  constructed  for  Meyer's  work  [3].  An  im- 
proved procedure  was  established  using  a  cylinder  constructed 
from  carbon  impregnated  paper.   Final  data  were  collected  on 
a  redesigned  paper  cylinder. 
1.   Meyer's  Cylinder 

The  following  is  a  brief  description  of  Meyer's 
cylinder.   A  detailed  description  is  contained  in  reference 
3.   The  cylinder  is  shown  in  Figure  1.   It  was  constructed 
of  acrylic  tubing  with  an  outer  diameter  of  4.47  inches. 
The  cylinder  was  32  inches  long  and  spanned  the  wind  tunnel 
test  section  from  roof  to  floor. 

A  surface  heat  flux  was  generated  by  energizing  the 
0.375  inch  wide,  0.003  inch  thick,  Nichrome  ribbon  which  was 
helically  wrapped  over  the  middle  15  1/8  inches.   Guard 
heater  circuits  prevented  end  losses  and  a  Silastic,  RTV 
foam,  limited  internal  convective  losses. 

Temperature  information  was  obtained  using  thermo- 
couples welded  to  the  inner  surface  of  the  Nichrome  ribbon. 
Two  pressure  taps  were  provided  for  obtaining  surface  pres- 
sure distributions.   The  cylinder  was  mounted  on  a  turntable 
that  could  be  rotated  from  outside  the  wind  tunnel.   This 
allowed  temperature  and  pressure  information  to  be  obtained 
at  all  angular  locations. 

In  the  present  set  of  experiments  the  cylinder  was 
installed  in  the  low  speed  Aerolab  wind  tunnel,  located  in 
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Figure  1  Meyer's  Test  Cylinder 
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the  Aeronautics  Laboratory  in  Halligan  Hall,  U.S.  Naval 
Postgraduate  School.   A  schematic  diagram  of  the  wind  tunnel 
is  provided  in  Figure  2. 

The  turbulence  level  of  the  wind  tunnel  was  measured 
by  Meyer  using  a  hot  wire  anemometer.   The  turbulence  level 
ranged  from  approximately  0.5%  to  0.7%  for  the  range  of 
Reynolds  numbers  studied  in  this  experiment. 

The  wind  tunnel  is  powered  by  a  lOOhp  electric  motor, 
and  has  a  four  speed  transmission.   The  tunnel  is  capable  of 
achieving  a  maximum  speed  of  200  mph  with  a  clear  test  sec- 
tion. 

In  order  to  obtain  temperature  information  on  Meyer's 
cylinder  using  liquid  crystal,  the  crystals  were  first  applied 
to  tapes  and  the  tapes  were  then  applied  to  the  cylinder 
surface  as  shown  in  Figure  3. 

The  tapes  were  made  by  spraying  Testors  flat  black 
paint  onto  Scotch  double  coated  tape  number  404.  A  diluted 
liquid  crystal  solution  was  then  applied  with  a  paint  brush 
in  two  coats.  Finally,  the  entire  assembly  was  sealed  with 
Polyurethane.  The  process  produced  tapes  with  an  approximate 
thickness  of  0.009  in.  and  a  crystal  layer  of  approximately 
.004  inches. 

As  noted  by  Meyer,  the  surface  of  this  cylinder  had 
numerous  surface  irregularities   caused  by  the  difference  in 
coefficients  of  thermal  expansion  between  the  acrylic  tubing 
and  the  Nichrome  ribbon.   These  irregularities  can  be  seen 
in  Figure  3,  and  were  the  primary  reason  for  designing  a 
new  cylinder  for  continuing  the  investigation. 
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Figure  3.   Test  Cylinder  with  Surface  Irregularities 
and  the  Liquid  Crystal  Tapes. 
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2.   Preliminary  Carbon  Impregnated  Paper  Cylinder  Design 

A  cylinder  that  did  not  exhibit  large  surface  irreg- 
ularities such  as  possessed  by  Meyer's  cylinder  was  construct- 
ed from  a  commercially  produced  resistive  paper  known  as 
Armstrong  Temsheet.   Temsheet  is  a  thin,  highly  flexible, 
electrically  resistive  material.   It  is  a  carbon  impregnated 
paper  containing  no  wires  or  ribbons.   The  nominal  thickness 
is  0.039  inches  and  the  electrical  resistance  is  approximately 
25  ohms  per  square.   The  heat  that  is  generated  when  a  con- 
stant current  is  passed  through  the  paper  is  uniform  to  with- 
in two  percent  from  point  to  point  over  large  areas. 

A  hollow  cylinder  with  an  outer  diameter  of  3.98 
inches  and  a  length  of  15  inches  was  formed  from  a  section 
of  Temsheet. 

To  hold  this  cylinder  in  place  two  10  inch  long,  4.5 
inch  outer  diameter  acrylic  tubes  were  attached  to  the  floor 
and  ceiling  of  the  32  inch  high  wind  tunnel  test  section. 
When  placed  inside  these  tubes,  the  Temsheet  cylinder  rested 
on  a  lip  inside  the  lower  cylinder.   Two  3.9  inch  diameter 
plugs  threaded  to  a  16  inch  center  rod  were  used  to  force 
the  Temsheet  snug  against  the  tubing  walls.   The  longitudi- 
nal seam  of  the  Temsheet  cylinder  was  sealed  with  tape  and 
was  placed  at  135°  from  forward  stagnation  on  the  side  not 
viewed  through  the  wind  tunnel  window.   The  details  of  con- 
struction can  be  seen  in  Figure  4. 

When  assembled,  a  12  inch  long  section  of  the  Tem- 
sheet cylinder  was  located  in  the  center  of  the  wind  tunnel 
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Figure  4.   Temsheet  Cylinder  with  Acrylic 
Tube  Holders. 
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test  section.   A  heat  flux  was  generated  by  using  a  Lamba 
Regulated  power  supply  model  LK345A  to  pass  current  through 
the  Temsheet. 

To  insure  uniform  application  of  voltage,  aluminum 
tape  was  used  as  electrodes  on  the  inner  surface  of  the  Tem- 
sheet.  The  tape  was  placed  circumferentially  on  the  surface 
2  inches  from  the  top  and  bottom  of  the  12  inch  test  section. 
Intimate  electrical  contact  between  the  tape  and  the  Temsheet 
was  established  by  applying  a  silver  based  conducting  paint 
on  the  edges  of  the  tape.   This  procedure  resulted  in  a  uni- 
formly heated  8  inch  test  section. 

To  guard  against  internal  free  convection  losses  the 
cylinder  was  filled  with  glass  wool. 

Since  Temsheet  is  black,  the  liquid  crystals  were 
applied  directly  to  the  test  surface.   To  insure  that  edge 
effects  would  not  interfere  with  readings,  no  liquid  crystals 
were  placed  within  one  inch  of  the  electrodes.   See  Appendix 
B  for  discussion  of  end  losses.   The  following  liquid  crystals 
bands,  3/8  inch  wide  and  separated  by  1/8  inch,  were  painted 
circumferentially  on  the  surface  of  the  Temsheet  cylinder. 
Reading  in  order  from  top  to  bottom:  R-49,  S-45,  S-43,  R-41, 
S-40,  S-38,  R-37,  S-36,  S-34,  R-33,  S-32,  and  S-30. 

The  surface  of  the  cylinder  was  marked  every  five 
degrees  to  allow  one  to  accurately  locate  the  various  iso- 
therms displayed  by  the  liquid  crystals. 
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3.   Final  Temsheet  Cylinder  Design 

The  protruding  lips  on  the  acrylic  tubes  which  held 
the  preliminary  Temsheet  cylinder  produced  an  undesirable 
secondary  flow  in  the  direction  of  the  cylinder  axis.  This 
axial  velocity  component  caused  the  overall  flow  to  be  highly 
three  dimensional.  Since  only  two  dimensional  flow  was  de- 
sired, a  new  set  of  supports  were  designed  that  exhibited  a 
constant  outer  diameter  to  the  flow.  The  final  cylinder  is 
shown  in  Figures  5  and  6. 

The  upper  and  lower  supports  for  the  Temsheet  test 
section  were  3.96  inch  diameter  wooden  cylinders  permanently 
attached  to  the  floor  and  ceiling  of  the  wind  tunnel.   The 
lower  cylinder  had  a  one  inch  diameter  hole  drilled  through 
its  axis  to  allow  the  passage  of  electrical  leads. 

The  wooden  cylinders  were  each  12  inches  long  with 
the  last  2  inches  on  the  ends  turned  down  approximately  .04 
inches.   This  allowed  the  Temsheet  to  be  attached  to  the 
cylinders  with  double  backed  tape.   The  outer  diameter  of 
the  assembly  was  then  constant  over  the  entire  wind  tunnel 
height. 

The  electrodes  were  attached  as  before  except  that 
the  test  section  was  only  6  inches  long. 

The  following  liquid  crystals  were  painted  directly 
on  the  Temsheet  surface;  recording  from  top  to  bottom:  R-4  9, 
S-45,  S-43,  S-40,  S-38,  S-36,  S-34,  and  S-32. 

A  second  cylinder,  utilizing  these  same  bases,  was 
constructed  with  one  crystal,  S-4  3,  covering  the  entire  sur- 
face from  electrode  to  electrode. 
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Figure  5.   Wooden  Bases  for  Final 
Temsheet  Cylinder. 
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Figure  6.   Final  Temsheet  Cylinder 
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Finally,  both  cylinders  were  marked  every  five  degrees 
and  filled  with  glass  wool  as  before. 

IV.   EXPERIMENTAL  PROCEDURE 

This  portion  of  the  investigation  consisted  of  three  dis- 
tinct phases. 

Initially,  the  cylinder  used  in  Meyer's  investigation  [3] 
was  used  to  check  the  liquid  crystal  readings  against  those 
of  the  installed  thermocouples.   The  procedure  outlined  by 
Meyer  was  used  to  bring  the  cylinder  to  steady  state.   The 
liquid  crystal  temperature  display  was  recorded  and  thermo- 
couple readings  taken  at  the  forward  stagnation  point  and 
then  every  15  degrees  of  arc. 

Some  difference  was  expected  between  the  two  temperature 
distributions  obtained.   This  difference  was  due  to  the  loca- 
tions of  the  two  sensors.   The  thermocouples  were  located  on 
the  inside  of  the  Nichrome  ribbon  while  the  liquid  crystals 
read  temperature  on  the  surface  of  the  cylinder.   Appendix  C 
contains  an  analysis  of  the  drop  in  temperature  between  the 
thermocouple  and  liquid  crystal  surface. 

Several  items  were  noted  during  this  phase  of  experimen- 
tation.  First,  the  degree  of  agreement  between  thermocouple 
and  liquid  crystals  depended  on  noting  the  point  at  which  a 
given  crystal  formulation  began  its  color  transition.   This 
was  not  suprising  since  calibration  procedures  mark  only  the 
beginning  of  color  transition  and  a  given  color  may  be  visible 
over  fractions  or  multiple  degrees  Fahrenheit. 
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Secondly,  the  surface  irregularities  caused  the  creation 
of  local  hot  and  cold  spots.   This  could  easily  be  seen  while 
the  cylinder  was  warming  up  or  when  the  cylinder  was  energized 
with  no  air  flow.   Further,  the  flow  pattern  in  the  area  of 
the  irregularities  was  altered.   This  was  visible  as  small 
regions  of  pulsating  liquid  crystal  colors  directly  behind 
the  irregularity. 

Additionally,  the  number  of  surface  irregularities  in- 
creased as  the  number  of  heating  and  cooling  cycles  on  the 
cylinder  were  increased. 

It  required  approximately  12  hours  for  Meyer  to  collect 
a  complete  set  of  thermocouple  data.   Using  the  liquid  crys- 
tals this  time  could  be  reduced  to  approximately  two  to  four 
hours.   Even  this  was  considered  too  long  for  practical  in- 
vestigations . 

It  was  decided  to  construct  a  new  cylinder  which  would 
have  a  smooth  surface  and  allow  data  to  be  collected  within 
a  reasonable  time  period. 

Before  proceeding  with  a  description  of  the  experimental 
procedure  using  the  Temsheet  cylinder,  a  brief  description 
of  the  temperature  distribution  that  exists  on  a  uniformly 
heated  cylinder  placed  in  a  cross  flow  of  air  will  be  given. 
This  description  together  with  the  accompanying  Figures  7 
and  8  should  prove  helpful  in  interpreting  the  liquid  crystal 
results. 

The  flow  of  a  real  fluid  is  best  analyzed  by  the  changes 
occurring  within  the  boundary  layer  that  forms  near  the 
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Figure  7.   Schematic  of  a  Cylinder  in  Subcritical 
Flow  with  the  Resulting  Temperature  and 
Heat  Transfer  Distributions. 
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Figure  8.   Schematic  Diagram  of  a  Cylinder  in 
Critical  Flow  with  the  Resulting 
Temperature  and  Heat  Transfer  Distribution 
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surface.   These  boundary  layers  are  classified  as  either 
laminar  or  turbulent.   The  flow  within  a  laminar  layer  is 
characterized  by  fluid  particles  moving  in  a  smooth,  orderly 
fashion  and  in  layers.   A  turbulent  boundary  layer  is  charac- 
terized by  individual  fluid  packets  flowing  in  irregular 
patterns  within  the  boundary  layer. 

The  heat  transfer  coefficients  and  temperature  distribu- 
tions are  functions  of  both  the  thickness  of  the  boundary 
layer  and  the  type.   As  the  boundary  layer  increase  in  thick- 
ness with  angular  location,  the  fluid  velocity  is  retarded 
resulting  in  an  increased  thermal  resistance  to  heat  transfer. 
In  the  case  of  a  uniformly  heated  cylinder  the  surface  tem- 
perature must  increase  to  maintain  the  wall  heat  flux  at  a 
constant  value.   The  nature  of  turbulent  flow  is  conductive 
to  better  heat  transfer. 

Flow  past  a  cylinder  may  be  classified  as  either  sub- 
critical  or  critical. 

In  subcritical  flow  the  laminar  boundary  layer  grows 
from  a  minimum  at  the  forward  stagnation  point  to  a  maximum 
at  an  angular  location  of  80-85  degrees.  At  this  point  the 
kinetic  energy  of  the  fluid  is  attenuated  sufficiently  that 
the  adverse  pressure  gradient  present  on  the  surface  can  no 
longer  be  overcome.  The  laminar  boundary  layer  separates 
from  the  surface.  The  surface  is  then  scrubbed  by  the  tur- 
bulent actions  of  the  wake. 

Figure  7  depicts  a  typical  subcritical  flow  pattern  and 
the  resulting  trends  in  both  temperature  and  heat  transfer 
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coefficient  with  angular  location.   Note  that  the  maximum 
temperature  exists  at  the  separation  point.   This  corresponds 
to  the  point  of  maximum  thermal  resistance,  or  minimum  heat 
transfer  coefficient. 

Critical  flow  is  characterized  by  an  increasing  laminar 
boundary  layer  thickness  on  the  forward  portion  of  the  cylin- 
der just  as  in  the  subcritical  case.   In  the  vicinity  of 
80-85  degrees,  however,  the  flow  begins  a  transition  from 
laminar  to  turbulent  flow.   The  point  at  which  transition 
begins  is  known  as  the  laminar  separation  point.   The  laminar 
boundary  layer  separates  but  the  kinetic  energy  of  the  fluid 
is  high  enough  that  the  flow  reattaches  to  the  surface  at 
some  point  aft  of  the  laminar  separation  point.   A  turbulent 
boundary  layer  then  develops  and  finally  separates  from  the 
surface  in  the  region  of  110-130  degrees  (the  actual  point 
of  reattachment  and  ultimate  separation  is  a  function  of  the 
Reynolds  number) . 

Figure  8  depicts  the  critical  flow  pattern  and  the  re- 
sulting trends  in  both  temperature  and  heat  transfer  coeffi- 
cient with  angular  location. 

The  development  of  a  consistent  method  of  obtaining 
accurate  heat  transfer  data  with  liquid  crystals  was  the 
objective  of  the  next  phase  of  experimentation.   The  tempera- 
ture fields  that  existed  on  the  Temsheet  cylinder  held  by 
acrylic  tubes  were  observed  under  steady  state  conditions  in 
the  wind  tunnel.   The  best  procedure  for  obtaining  accurate 
information  was  developed  by  analyzing  the  visual  temperature 
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fields  and  comparing  them  with  those  of  the  previously  dis- 
cussed theory  of  angular  temperature  and  heat  transfer  dis- 
tributions . 

As  noted  earlier,  the  most  accurate  temperatures  were 
found  at  the  beginning  of  color  transitions  on  the  liquid 
crystals.   In  areas  of  large  temperature  gradients  there  was 
no  problem  in  determining  the  start  of  red,  green,  or  blue. 
The  particular  liquid  crystal  in  the  temperature  range  at 
these  points  would  undergo  the  red  to  green  to  blue  transi- 
tion in  a  short  interval.   The  point  of  blue  transition  was 
easily  noted. 

The  point  of  transition  in  the  area  of  stagnation  was 
much  more  difficult  to  observe.   The  temperature  gradients 
in  this  region  were  very  shallow.   It  was  not  uncommon  for 
one  liquid  crystal  to  exhibit  approximately  the  same  color 
over  15  to  20  degrees  of  arc.   The  actual  transition  to  this 
color  was  usually  not  apparent.   The  error  in  reading  a  tem- 
perature in  this  region  could  be  as  large  as  1°F.   To  obtain 
a  precise  temperature  the  beginning  of  a  transition  was 
forced  by  adjusting  the  power  supply  until  a  previously  un- 
changed liquid  crystal  just  began  its  transition.   This  was 
verified  by  decreasing  the  voltage  a  few  tenths  of  a  volt 
and  observing  the  liquid  crystal  return  to  its  previous 
color. 

If  a  clear  cut  transition  from  red  to  green,  or  green  to 
blue  occurred  in  any  region  of  shallow  temperature  gradients, 
it  was  considered  precise  and  used.   This  type  of  transition 
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usually  occurred  in  angular  regions  15  to  60  degrees  aft  of 
the  forward  stagnation  point. 

The  hottest  point  on  the  cylinder  surface  was  at  separa- 
tion (Figure  7) .   This  point  was  found  by  adjusting  voltage 
until  the  first  point  of  red  color  was  observed  on  an  un- 
changed crystal  or  decreasing  power  on  one  that  was  already 
through  its  transition  until  the  very  last  red  point  remained, 
The  location  was  checked  by  repeating  the  procedure  on  a 
different  liquid  crystal  band.   The  two  points  invariably 
fell  on  a  vertical  line. 

By  combining  the  procedures  repeatable  results  were  ob- 
tained for  a  given  Reynolds  number. 

While  experimenting  with  the  above  procedure  the  Reynolds 
number  was  maintained  at  approximately  50,000,  clearly  in  the 
subcritical  flow  regime.   However,  a  local  cold  spot  not  at- 
tributable to  normal  subcritical  flow  patterns  existed  in  the 
region  of  90  degrees.   Small  strings  placed  on  the  cylinder 
confirmed  that  undesirable  secondary  flow  patterns  were  gen- 
erated by  the  lips  on  the  acrylic  tubes.   Due  to  this  flow 
anomaly  a  third  phase  of  experimentation  employing  an  im- 
proved wooden  base  was  initiated.   The  objective  of  the  third 
phase  of  the  investigation  was  to  gather  data  for  comparison 
with  conventional  works.   This  was  done  in  three  steps.   The 
new  wooden  base  Temsheet  cylinder  was  tested  for  two  dimen- 
sional flow  characteristics  and  edge  losses.   The  roughness 
of  the  Temsheet  and  liquid  crystal  coatings  were  examined. 
Data  for  Reynolds  numbers  ranging  from  38,000  to  148,000  were 
obtained  for  analysis. 
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For  data  collection  the  Temsheet  cylinder  with  eight 
liquid  crystal  bands  was  installed  in  the  wind  tunnel. 

The  airstream  thermocouple  and  its  reference  ice  bath 
were  inspected  and  temperature  recorded  for  initial  calcula- 
tion of  airstream  properties  and  velocity. 

A  U-tube  manometer  and  micromanometer  used  for  airspeed 
indication  were  carefully  zeroed. 

The  power  supply  was  energized  and  voltage  adjusted  to 
cause  all  liquid  crystal  bands  to  undergo  a  transition. 
This  served  two  purposes.   First,  if  the  glass  wool  packing 
was  not  installed  snuggly  a  point  of  temperature  discontinu- 
ity would  appear.   The  packing  could  then  be  adjusted  prior 
to  establishing  air  flow.   Secondly,  the  cylinder  could  be 
preheated  to  reduce  the  time  necessary  to  reach  steady  state. 

The  wind  tunnel  was  started  and  the  air  speed  was  in- 
creased using  the  previously  calculated  pressure  drop  on  the 
U-tube  manometer. 

Since  the  U-tube  manometer  was  not  very  sensitive  over 
the  entire  range  of  Reynolds  numbers,  especially  at  the  lower 
Reynolds  numbers,  the  pressure  drop  read  on  the  micromanometer 
was  used  to  recalculate  air  speed  and  Reynolds  number.   This 
proved  to  be  a  critical  part  of  the  experiment. 

Power  was  increased  until  the  upper  liquid  crystal,  R-49, 
began  its  red  transition  at  separation.   Power  was  continu- 
ally adjusted  to  maintain  red  at  separation  until  steady 
state  was  reached.   This  usually  occurred  within  fifteen 
minutes  and  was  typified  by  a  period  of  approximately  five 
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or  more  minutes  of  a  constant  red  without  adjusting  power. 
This  was  checked  by  decreasing  power  a  few  tenths  of  a  volt 
and  noting  the  red  spot  disappear. 

Data  were  collected  after  restoring  power  to  its  previous 
setting.   The  angular  location  of  each  liquid  crystal  color 
transition,  voltage,  air stream  temperature,  and  micromanom- 
eter  height  were  recorded. 

For  critical  flows  the  procedure  for  locating  separation 
was  followed  to  locate  the  points  of  laminar  separation, 
reattachment,  and  final  separation. 

The  test  for  edge  effects  and  two  dimensional  flow  char- 
acteristics consisted  of  placing  a  Temsheet  cylinder  coated 
with  liquid  crystal  S-43  from  electrode  to  electrode  in  the 
wind  tunnel. 

The  wind  tunnel  was  started  and  the  cylinder  heated  to 
steady  state.   The  isotherms  were  then  observed  for  straight- 
ness  and  end  effects. 

Although  the  surface  of  the  Temsheet  was  smooth  to  the 
touch  the  actual  surface  roughness  and  effects  of  the  liquid 
crystals  were  not  known.   A  piece  of  Temsheet  coated  with 
liquid  crystals  was  examined  under  an  electron  microscope. 
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V.   RESULTS 

A.   TEMPERATURE  DISTRIBUTION  AND  HEAT  TRANSFER 

The  numerical  data  and  data  reduction  technique  are  con- 
tained in  Appendix  A.   Briefly,  the  local  heat  transfer  co- 
efficient was  obtained  by  dividing  the  heat  generated  per 
unit  area  by  the  temperature  difference  between  the  surface 
and  airstream  temperatures. 

h  =  q/AT 
The  value  of  q  was  determined  by  dividing  electrical  power 
generated  by  the  surface  area 

q  =  (V2/R) 
A 

The  local  convective  heat  transfer  coefficient,  h  ,  was  then 

c 

obtained  by  subtracting  the  radiation  coefficient,  h  ,  (see 
Appendix  A)  from  h. 

hc  =  h-hr 
The  Nusselt  number  (Nu)  was  calculated  using  the  expres- 
sion 

Nu  =  hcD 

k 

and  the  Froessling  number  calculated  from 

Fr  =  Nu 


Data  are  presented  graphically  in  this  section.  Figures 
of  temperature  distribution,  Nusselt  numbers,  and  Froessling 
numbers  are  included. 
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Figure  9  shows  a  sketch  of  a  typical  subcritical  flow 
temperature  distribution  and  a  representation  of  the  isotherms 
seen  on  the  liquid  crystal  bands  on  the  cylinder  surface. 

Critical  flow  was  noted  at  a  Reynolds  number  of  approxi- 
mately 120,000.   The  shape  of  a  typical  critical  flow  temper- 
ature distribution  was  shown  in  Figure  8.   Figure  10  is  a 
color  photo  of  the  actual  temperature  distribution  on  the 
cylinder  surface  for  this  critical  flow  condition.   Note  that 
on  the  upper  liquid  crystal  band,  R-49,  the  presence  of  the 
separation  region  is  characterized  by  a  maximum  temperature 
at  115°  (denoted  by  a  red  patch)  with  lower  temperatures  on 
either  side  (denoted  by  black) .   The  blue  region  forward  of 
this  region  indicate  higher  temperatures  and  bracket  the 
laminar  separation  point. 

The  Nusselt  numbers  for  several  Reynolds  numbers  are 
plotted  in  Figure  11.   Several  items  should  be  noted  from 
these  curves. 

First,  the  general  shape  of  the  curves  is  consistent 
with  the  previous  work  of  Giedt  and  Meyer  [2,3]  in  both  crit- 
ical and  subcritical  flows. 

The  points  of  minimum  heat  transfer  coefficients  occur 
in  the  region  of  8  0  to  90°.   However,  unlike  the  work  of 
Giedt,  these  points  move  forward  on  the  cylinder  with  in- 
creasing Reynolds  numbers.   This  may  be  due  to  a  slight 
deformation  of  the  cylinder  as  velocity  increases.   Deforma- 
tion was  not  visible  to  the  naked  eye  at  Reynolds  numbers 
below  150,000.   It  was  interesting  to  note  that  the  isotherms 
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Figure  9.   Liquid  Crystal  Temperature  Field  and  Resulting 

Temperature  Distribution  on  the  Cylinder  Surface 
Subcritical  flow  conditions  indicated. 
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Figure  10.   Liquid  Crystal  Temperature 
Field  in  Critical  Flow. 
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Figure  11.   Heat  Transfer  Results  at  Reynolds  Numbers 
of  38,000,  88,000,  and  148,000. 
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depicted  on  the  cylinder  surface  were  no  longer  vertical 
when  deformation  was  visible.   This  can  be  seen  in  Figure  12. 

The  Froessling  numbers  for  the  laminar  flow  regions  are 
plotted  in  Figure  13  for  comparison  with  the  theory  of  Schuh 
[1] .   Comparison  is  within  experimental  uncertainty  near 
stagnation  and  then  drops  off  near  separation.   The  latter 
trend  is  consistent  with  the  work  of  Giedt  and  is  most  prob- 
ably explained  by  the  fact  that  Schuh  used  an  ideal  pressure 
distribution  as  opposed  to  the  distribution  encountered  in 
viscous  flows. 

B.  EDGE  EFFECTS 

Figure  14  showed  the  isotherms  on  the  cylinder  coated 
with  one  liquid  crystal,  S-43.   They  clearly  indicate  that 
heat  is  lost  near  the  top  and  bottom  edges  of  the  coated 
region. 

The  edge  effects  are  indicated  by  a  slight  curvature  of 
the  isotherms  near  the  electrodes.   This  is  expected  from 
the  model  developed  in  Appendix  B  and  was  the  reason  for  not 
taking  data  within  an  inch  of  the  electrodes. 

C.  SURFACE  ROUGHNESS  OF  TEMSHEET 

Figures  15  and  16  show  the  Temsheet  liquid  crystal  inter- 
face magnified  121  and  27  times  respectively.   The  point  to 
note  is  that  the  liquid  crystals  fill-in  irregularities  in 
the  surface  and  appear  to  be  20-30  micron  spheres  as  speci- 
fied by  the  manufacturer. 
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Figure  12.   Cylinder  Deformation  at  Reynolds 
Numbers  Greater  than  150,000. 
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Figure  13.   Comparison  of  Experimental  Results  of  the 
Present  Investigation  with  the  Theoretical 
Solution  of  Schuh. 
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Figure  14.   Isotherms  Showing  End  Losses  and 
Two  Dimensional  Flow  Pattern. 
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Figure  15. 


Temsheet-Liquid  Crystal  Surface 
at  121X. 


Figure  16. 


Temsheet-Liquid  Crystal  Surface 
at  2  7X. 
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Fage  and  Warsap  in  reference  18  predict  a  transition  to 
critical  flow  at  a  Reynolds  number  of  approximately  148,000 
for  flow  over  a  cylinder  with  a  surface  roughness  equivalent 
to  the  roughness  of  the  crystal  coated  Temsheet  (k/D=2xlO-1*)  ; 
k  being  the  critical  roughness  height  and  D  the  diameter  of 
the  cylinder.   The  actual  transition  occurred  at  approxi- 
mately 120,000. 

VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  liquid  crystal  thermographic  technique  developed  in 
this  investigation  provides  an  excellent  means  of  obtaining 
both  qualitative  and  quantitative  heat  transfer  information 
on  heated  objects  placed  in  forced  convection  environments. 
Using  the  technique  it  was  possible  to  quickly  and  easily 
obtain  information  on  the  variation  of  the  Nusselt  number 
around  the  circumference  of  a  uniformly  heated  right  circular 
cylinder  placed  in  a  cross  flow  of  air.   The  technique  also 
allowed  one  to  visually  observe  the  effects  of  flow  separa- 
tion, the  turbulent  boundary  layer,  and  the  turbulent  wake 
on  the  surface  temperature  of  the  cylinder.   Colored  movies 
taken  of  a  cylinder  coated  with  a  single  liquid  crystal  are 
especially  vivid  in  their  display  of  the  influence  of  the 
turbulent  wake  on  the  cylinder  surface  temperature.   In  the 
wake  region  the  crystals  alternately  dim  and  glow  in  response 
to  the  "scrubbing"  action  caused  by  random  bursts  of  cool 
fluid  impacting  on  the  surface. 
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The  ability  to  visually  observe  turbulent  flow  pattern 
effects  on  surface  temperatures  suggests  an  excellent  means 
for  studying  the  influence  of  free  stream  turbulence  on  the 
heat  transfer  rates  of  heated  objects.   However,  in  order 
that  such  a  study  be  conducted  in  a  quantitative  manner,  the 
thermal  response  time  of  the  liquid  crystals  must  first  be 
determined.   Parker  [19] ,  using  a  capacitor  discharge  tech- 
nique, found  that  a  0.001  inch  film  of  liquid  crystals 
coated  on  a  thin  stainless  steel  foil  responded  in  approxi- 
mately 0.036  seconds  to  a  step  change  in  foil  temperature. 
It  is  possible  that  a  thin  coat  of  liquid  crystals  placed  on 
a  material  such  as  Temsheet  would  respond  even  faster  than 
this  due  to  the  fact  that  the  crystals  are  partially  ab- 
sorbed by  the  Temsheet.   The  ideal  situation  would  be  if  the 
liquid  crystals  responded  at  the  same  rate  as  the  Temsheet 
itself.   An  experiment,  similar  to  the  one  conducted  by  Parker, 
needs  to  be  conducted  to  determine  response  time  information 
on  the  liquid  crystals. 

As  a  final  recommendation,  the  phenomenon  shown  in  Figure 
17  is  offered  as  a  possible  topic  for  future  investigation. 
During  the  final  phase  of  collecting  wind  tunnel  data  for 
the  present  investigation,  it  was  noted  that  the  cylinder 
coated  with  liquid  crystal  S-43  displayed  alternate  hot  and 
cold  spots  along  the  separation  line.   These  spots  were  uni- 
formly spaced  and  seemed  to  be  caused  by  some  flow  phenomenon, 
perhaps  a  series  of  vortices  such  as  one  observes  in  flow  in 
a  curved  channel.   Whatever  the  cause,  the  hot  and  cold  spots 
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Figure  17.   Hot  and  Cold  Spots  on  the  Cylinder 
Surface  along  Separation  Line. 
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definitely  existed  as  is  readly  apparent  in  Figure  17.   Pre- 
cise measurements  were  not  taken  at  the  time  the  phenomenon 
was  observed  (the  Reynolds  number  was  approximately  75,000) 
and  no  explanation  is  offered  for  its  existence.   It  should 
be  noted,  however,  that  such  a  phenomenon  may  well  have  gone 
undetected  if  thermocouples  were  used  as  the  temperature 
sensor.   Additional  research,  perhaps  using  liquid  crystals 
for  temperature  sensing  and  smoke  for  flow  visualization, 
needs  to  be  conducted  to  explain  the  observed  phenomenon. 


47 


ft- 


■f    '.-'. 


3B3i 

b 


.  _  - 
......       • 


OS" 


APPENDIX  A 
DATA  AND  DATA  REDUCTION 

In  order  to  experimentally  determine  the  Nusselt  number 
as  a  function  of  angular  location  on  the  right  circular 
cylinder  under  investigation,  a  relationship  between  the  con- 
vective  heat  transfer  coefficient,  hc ,  and  the  surface  heat 
flux  produced  by  the  Joulean  heating  effect  in  the  Temsheet, 
V2/RA,  was  needed.   This  relationship  is  easily  developed  by 
performing  a  simple  energy  balance  on  an  elemental  volume, 
see  Giedt  [2]  or  Meyer  [3] .   The  result  is 

V2   kt  dT2 

hc  -  RA  +  rQ2  d92  -  hr 

(T-TJ 

where 

h„  =  Convective  heat  transfer  coefficient 

V  =  Voltage  impressed  across  the  test  section 

R  =  Electrical  resistance  of  test  section 

A  =  Area  of  test  section 

t  =  Thickness  of  test  section 

k  =  Thermal  conductivity  of  test  section  material 

r  =  Radius  of  the  cylinder 

T  =  Temperature  at  angular  location  9 

0   =  Angular  location 

T^  =  Air  stream  temperature 

hr   =   Radiation  heat  transfer  coefficient  -  aF._,{T+T  }{T2+T  2} 

IT  \      «       oo  00 

Fj_2  =  Radiation  exchange  factor  between  cylinder  and  sur- 
roundings 
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a     -   Stefan-Boltzman  Constant 

An  exact  value  for  thermal  conductivity  of  the  Temsheet 
material  was  not  available,  but  assuming  that  its  value  is 
approximately  the  same  as  other  carbon  paper  products,  the 
product  kt  is  calculated  to  be  .000195  BTU/hr-°F.   Since  this 
value  is  so  low  it  was  felt  that  the  conduction  term  in  the 
expression  for  convective  coefficient  could  be  neglected. 
The  expression  used  for  the  convective  coefficient  then  takes 
the  final  form: 

hc  =  Yl "  hr 

RACT-TJ 

The  Nusselt  number  is: 

Nu  =  2hcro 
kair 

The  Froessling  number  is: 

Fr  =  Nu/  T|Re" 
where 

Re  =  Reynolds  number  =  V^D/v 
Voo  =  Air  velocity 
v  =  Air  kinematic  viscosity 

A  sample  calculation  is  provided  for  illustration  using 
the  following  values  from  the  data  at  a  Reynolds  number 
47,818: 
9   =  60* 
e   =  .9 

a      =  .1714  x  10~8BTU/hr-ft2-oR" 
^  =  62.7°F 
T   =  96.3°F 
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Voo  =   22.5   ft/sec 
V     =23.7  Volts 
R     =  12.0   a 

D/K  =  22  hr-'F-ft» 
BTU 

Re  =  47,818 

h  = (23.7V)2 

(12.0ft)  (.52ft*)  (33.6T)  ( .  293WATTS/BTU/hr 

h   ■  9.2  BTU 

hr-ft2-°F 

h  -  .9(.1714  x  lO-8BTU/hr-ft2-°R*)  (556.3+522.7) 
(556.32f552.72) °R3 

hr  =  .97  BTU 

hr-ft2-°F 

h   =  8.2  BTU 

hr-ft2-°F 

Nu  =  (8.2  BTU       )  (22  hr-°F-ft2) 
hr-ft2-°F     BTU 

Nu  =  180.4 

Fr  =  180.4     =  .82 
^47,818 

The  data  sheets  for  all  seven  Reynolds  numbers  investi- 
gated follow. 
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DATA 


A1r  Temp.  =  6£>~°F;        /    =  2°f°  f?£ 

v  ft 

Air  Speed  =    -/?  cmH20 
■   ft    ft/sec 

Re    =     365oo 

Correction  Factor  =  /-02&7 
Recorr=    37,*e.T" 
Voltage  =  3*6" 4  V 
Heat  Flux  =  ^r?  BTU/hr.ft2 


Date:  /-**''?¥ 
Time:       /£,</£~~ 
Crystal  reader:  T^/ui 
Data  sheet:  (2c0/o£^ 
Cylinder  Dia.  =  '33>  /V 
Cylinder  Area  =  -  S~>  / 
Resistance  =  '^  ~^— 

Crystals:  ^  „.  _ 

^^        s  3^ 

5-/^     ssy 


Film  Temp  ■    ?£>  °F; 

0/     _^/^hr-°F-ft2 

S</0 

53  a 

8 

CRY. 

TEMP 

r°Fi 

AT 

h 

hr 

he 

Nu 

Fr 

Comments 

0 

5  3S-G- 

/02.<l 

3*.* 

9-73 

Aoo 

g.73 

/76 

.f/ 

-     /f->y.0 
^  "5 £/}#/£#/  /on/ 

60 

5V6'8 

/07,/ 

7/- 6 

SmS* 

/.Of 

7,C( 

''6& 

«J5 

71 

S938 

//a.  3 

76,£ 

Z67 

7^7 

6^</ 

/VS~ 

0  75~ 

7V 

59tT/S 

//£~,/ 

776 

7=23 

7c<y 

&,/<? 

/Ss~ 

.^7 

77 

/?v^ 

W>7 

££^a 

6.3^ 

/^<j<r- 

5-  3  y. 

//£ 

a    6 

/"3 

/?#/3 

/7/-7 

5~£,a 

6.57 

^7  <r 

<5~'3tj 

//6 

'^ 

/30 

5</5-a 

//5~»/ 

77  £ 

"7^3 

/.CU 

C/9 

/35- 

./37 

5Y3/3 

/y«?«3 

76,7 

^47 

76  3 

Ca    -&Q 

/75~ 

/CAT 

S</£>6 

/07,/ 

77  4 

5-6^ 

AC 7 

7.6/ 

/££ 

/S3- 

53gV3 

/a&,9 

2?V 

&5;F 

Aoq 

ST'Sj- 

/E^~ 

/*.r- 

6^:  * 

994 

53.  ^ 

/asr 

'2F 

9* 60 

^/9 

/So 

S3V6- 

?SZJt 

59,7 

-£?#,? 

^7 

//  // 

3</^ 

e/// 

1/U6/£Jb 

/'ex. 

~/J<5^- 

.     t/  = 

,5P,/ 

\?        /& 

'  =  65~6 

5-6 

53  7 /3 

9&,~2 

3o.% 

5.  £7 

,?i 

C~7.G? 

-*<S 

c$- 

S3£/^ 

99./ 

S3,9 

7.8$ 

<?7~ 

G'S9~ 

/S~0 

-77 

7S- 

S3?/s 

/T5.9 

77.  V 

6,7  5- 

/■Co 

-7T55- 

'<#,7 

*£7 

7% 

S?2>/3 

/67,/ 

V/C 

6^77 

Ac/ 

67// 

s/fi 

.<^/ 

#3. 

5^3,^ 

//a.  3 

#64 

5:?/ 

/.0  3 

7-6  7 

/03. 

.^~^ 

s-tf 

«stf-r<e 

7/5"  / 

77  6 

S-,35- 

/,<j</ 

/'3  7 

9s~ 

-77 

5>7 

/*?#/? 

/£?// 

7<6^ 

#75- 

/,a<T 

3,7- 

Sr/ 

^7A 

/<0  0 

5YiTyfi 

//y./ 

99 -c 

S73£ 

/■oq 

^•3<7 

75' 

/o^- 

5V-?/* 

S/3.S 

</£.$ 

ST  7/ 

4>>> 

X6^ 

//X 

/M 

5^5/5 

/C7,/ 

9/.C 

6,</,5 

/.c/ 

677/ 

//$ 

/3H 

S?.?£ 

/03/7 

5:<,v 

«^<r 

/<or 

5^7  <- 

/Sn 

/v<c 

5  3fA 

??,« 

'•7  7 

7^f 

■99 

<r.  56 

'  s~o 

rS-<r 

53V/? 

9/,S 

3*,£ 

■-^7 

<9J 

Z<;9 

/6§> 

/£<r 

SSifr 

67  6 

7/./ 

7.77- 

'9/ 

S-Sfi 

7>  / 

1 
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DATA 


AlrTemp.  *&?.7  °F;      D/y  =  2070 
Air  Speed  =  6V23  anH20 

Re    =    AU,S7S- 
Correction  Factor  =  7,02br] 

Re  =    47,?/? 

corr  ; 

Voltage  =  23.7  v 

Heat  Flux  =  308  BTU/hr#ft2 


sec 
ft 


Date:    /-24-7? 
Time:     703  OO 
Crystal  reader:   FIELD 
Data  sheet:      Coope/2 
Cylinder  Di a.  =   O.  33c?  fir 
Cylinder  Area  =  O.  S/^fiT2- 
Resi  stance  =    J2.QJI 
Crystals: 


Film  Temp  =  Q5~  °F; 

D/     _??  hr-°F-ft2 

*4o 

S3* 

8 

CRY. 

TEMP 

r°Fi 

AT 

h 

hr 

he 

Nu 

FR 

Consents 

60 

S-248 

36.3 

33^6: 

$J  7 

^7 

3.1- 

/£0.  v 

.$2 

su/sc/irtc/lL  Eic  > 

NO    /30&ZLE 

72 

S  -ld>  6 

3-3.4 

36.7 

S.sD 

.?$ 

7,9/ 

J/3 

,74s 

76 

S3§_8Ud3.3 

41.2 

749 

-y? 

C,i/9 

/</1 

a    6  5" 

73 

S  40  B 

107.  / 

44  A 

(o.34 

/■DO 

ft* 

73/ 

o     Co 

$2 

34  2  B 

J/2. 3 

43.  & 

{..21 

70/ 

f.% 

//Q 

-52 

94 

s4£& 

/JS~J 

52  4 

SS8 

/■or 

V,ft£ 

707 

.*? 

27 

J143R 

/ '2/  7 

53.  O 

S.22 

/,o4 

</./y 

92- 

.   37 

^S'B 

7/5,  1 

S2.4 

f.88 

J.C7- 

t/.te 

707 

102 

5438 

J/ 2. 3 

43.  u 

(a-21 

J  .01 

S'.U- 

//y 

J/S 

S40/3 

107.1 

44.4 

(,.34 

/.CO 

f.f* 

/3/ 

/Z2 

S28S 

703.3 

412 

748 

.49 

ci-4 

/</3 

732 

S3GS 

33.4 

3L7 

£.33 

.4$ 

*7,¥/ 

/O,  3 

/40 

534  B 

3&.3 

33.  {a 

3/7 

■77 

%-* 

/$0.tj 

/5V 

SJZg 

32. 6 

23.3 

/Q.Sd 

•*(, 

7-l-T- 

2075- 

€//, 

9A/CVT/)   V 

'OirAGEt 

'0   27  6 

V  .  o  = 

m 

r^-- 

-62.7 

0 

538G 

102.  4 

30.7 

/o.so 

-.  p 

%s-i 

2C9 

-95 

u 

S  40  8 

/O  7. 1 

444 

9,  ^3 

/■CO 

4-  if 

7SS 

°$Y 

n 

S41& 

//a -3 

49.£> 

%4l 

/.oj 

7-<7 

/63 

o7<7 

74 

545 B 

/itS- 1 

5Z-4 

7.  90, 

/.ol 

6.9* 

7S3 

c  7 

77 

R10S 

1^.7. 

c/.s- 

0-79 

/<C± 

r-73 

22  £ 

<5% 

104 

/W/3 

/7Y-2 

6/-  S~ 

C-  ■  72 

/.of 

s:/3 

/20> 

120 

$4$$ 

//S./ 

&3.V 

7.9  4 

/-0% 

£-^ 

/-S3 

/3i 

543  3 

//■J-  3 

H9-L 

-7-7/ 

/iO/ 

7-70 

/C-3 

742 

S10S 

/o7/ 

44.4 

c/-  ?7 

/ .  06 

4^9 

/-33~ 

/SO 

S38& 

/Si,  </ 

4/. -p. 

/C./2- 

■99 

9/3 

20/ 

/S7 

3  563 

?£V 

3C  7 

// 3  4 

■797 

/S3t 

23$ 

/SO 

5  34G 

<75l% 

39  ■  S~ 

/  2  .  S  3 

-7<c 

/AS'/ 

'V,  / 

so0 

5  38B 

io'S3 

41. 1 

JO.  12 

c/9 

/.  /3 

2,0/ 

»77L 

52 


DATA 


A1r  Temp.  =  <cM  °F; 

Air  Speed  ■  .  C    cmfoO 
■  S3,  ft/sec 
Re    =  Cc,j2<fo 
Correction  Factor  =   7-03-&7 

Re  =  &S,oo2 

corr 

Voltage  =  2S-*r  V 

Heat  Flux  =  ^/7BTU/hr.ft2 


/    =2070 
v 


sec 
ft 


Date:     /-21-7V 
Time:     //y3~ 
Crystal  reader:  /r/^/7> 
Data  sheet:   Coop** 
Cylinder  Dia.  =  .  33  $t 
Cylinder  Area  =  <>s~/9  f'-A 
Resistance  =  /s.o-sl. 
Crystals: 


6Q2 


S3* 

say 


Film  Temp  = 

??°F; 

D/     = 

7/^hr-°F-ft2 
z/.x     RT|| 

!5</0 

5"S3- 

e 

CRY. 

TEMP 
(°F) 

AT 

h 

hr 

he 

Nu 

FR 

Conments 

0 

S^>%G 

/62,<t 

3f -O 

/^-33 

-9f 

/A  29 

2-</5~ 

<9u 

<-  CMAAJ6£J> 
yct-TAGS 

Too  -OZ-V0/- 

ItT 

538  B 

lO'icl 

<#5.  S- 

//•  75- 

<9? 

70.79 

3-  3s" 

»<? 

6"/ 

SVaG 

Jrr.  <l 

#?<7 

//■/7 

/.CO 

70,77 

^3  2 

*$$- 

sy 

s~&/3 

to) ,/ 

V^-7 

/OTJSL 

/■oo 

9^72. 

<5/<J. 

°%* 

?o 

SViS 

//7-3 

£#<£ 

9-?.r 

/■o^ 

S.7^ 

/9cJ 

»73 

y^ 

Svs~3 

/7&~,  / 

5-/  7 

9-5,  3 

/■0<7 

?</? 

'79 

■  09 

77 

/?^73 

/$</.& 

Q:-,<7 

77^- 

7-ocr 

6,9 

7<Jf 

-5-7 

./O-z, 

/ev^/3 

/2</-£ 

6>or% 

?,*s~ 

/■o$- 

^9 

/c/79 

7^0 

5V3-/5 

//5'/ 

6T7.  7 

7;i3 

y-0(7 

<?//T 

779 

/2  6 

5/^/5 

7S3,2, 

s7S<5 

9-7<r 

AOX 

S-77 

7?d 

/3<r 

5^/3 

/C?,/ 

(73.7 

7C$*L 

7.0 0 

992- 

£/<9 

/<7<r 

538-^6 

/ci-9 

</0,<r 

'7<7Z 

<99 

70.7? 

7  3  5- 

fsr*r 

S3G  /3 

997 

Be,.® 

03^<r 

,95 

7£<37 

J?^7 

/£5T 

^3y/? 

?£,3 

32,9 

'ts- 

'97 

73.$ 

2oy 

/So 

S  3^-G- 

93-6 

Zf.3. 

7£,3<7 

,  96 

/S~*3/f 

?3J- 

<r$- 

S3-/S 

96.^ 

'55.  f 

/ O.ty 

■  97 

9  7 

Z// 

<S7 

67 

337  A 

99-<A 

3<£ 

9.7^- 

■  99 

S-7? 

'?/ 

-  7^ 

70 

SaS^ 

m.5 

^5- 

S,i7 

■99 

7-&X 

7S7 

*  &Y 

/V 

s%>/? 

7c  7/ 

V^7 

&<o^ 

/•G<\ 

>,*R 

7S~3 

cS9 

76 

84^3/8 

/^■3 

V77 

y./# 

7<o% 

£'/& 

/3<7 

c  5 1 

7% 

S^5"yS 

7/S~,t 

6-77 

£>.79 

/■OU 

$~7^- 

/75~ 

*£ 

^5 

>fV?/-? 

727.7 

5~V?«3 

0-07 

7o<s 

^■97 

S69 

/CO 

5^5"/Q 

//S"// 

57-7 

£,7^ 

&</ 

S~7b" 

775- 

70  7- 

3V8t& 

vza 

x77  7 

7.S8 

/.07- 

6  ,/C- 

/3<7 

/OS 

S7079 

/07,7 

V7/ 

8.0& 

7-  Od 

ZC7 

'SX3 

7J*7 

S  36  & 

99(7 

.7^ 

7/r 

'99 

9  77 

79  / 

/3o 

S$i/8 

<?7.3> 

i'j?,  7 

/d<67 

■97 

9.7 

7// 

/</.? 

5  3*  A 

?2,& 

,*M 

.-7  OA 

-'/& 

//.  0£ 

9/77 

53 


DATA 


A1r  Temp.  =  ci.l  °F;      D/    -   loco  ||c 

v  ft 

Air  Speed  =  ' '**  cmH20 
■  &i   ft/sec 

Correction  Factor  ■  />oz&l 

Re„„„  =  A??.  ^65" 
corr      ' 

Voltage  ■  59-5  V 

Heat  Flux  =  57?  BTU/hr#ft2 


Date:  i'tt^Tt. 
Time:  /3as" 
Crystal  reader:  /-~^<-*> 
Data  sheet:  G-oop&z 
Cylinder  Dia.  ■  .  33.  Pt 
Cylinder  Area  =  ,<sv9#4 
Resistance  =  /**o 
Crystals: 

S<73  S3</ 


Film  Temp  =    S^  °F; 

0/    -2/.J?hr-»F-ft2 

5^c 

S3&. 

e 

CRY. 

TEMP 

AT 

h 

hr 

he 

Nu 

FR 

Conments 

0 

S3SG- 

/<?2,y 

38 /i 

/s~,/z. 

*f* 

/Y/  /S 

3C7 

c^y 

{,/a.c-r/PG-^ 
!Z>  -0,7/ 

?  -  «*? 

A/0   /204£6£ 

a£ 

532  8 

JC&.q 

3^.5 

Mrr 

.^f 

/3,3~6 

2L9C 

.96 

6~6> 

^/og- 

/££./ 

42.C 

/a  79 

//00 

/H.  79 

2  79 

.  %S 

60 

5^0(3, 

/o7.i 

V3 

/3,</  7 

/•  CO 

/2>Q7 

27  z 

'S3 

7o 

SU3B 

//S/3 

VS.5- 

/3<ai 

/■6/ 

//■Oo 

Ai/o 

.73 

V3 

^¥<6 

//SW 

57.0 

/ABS 

//02. 

/c,  33 

a^.y- 

,0,9 

?7 

f?*fl/& 

/7/.p- 

£4/ 

?<£> 

/C5' 

y-fr 

/9G 

'S-7 

/02- 

/979s 

/5^.ja 

<&?,/ 

9-C 

A'OfT 

^S-5- 

/?£ 

/0<T 

5YS-/S 

//<,-/ 

-57'CJ 

//.*<r 

AOJ- 

/t?.  5  5 

^s<r 

7/7- 

SifeB 

//2'3 

VZ-i- 

72 -6/ 

Ao/ 

//■  00 

TiMo 

/&2- 

J^C£ 

/07,/ 

^3 

/3,</7 

/-on 

7^7? 

277 

«6— 

S88/3, 

#8,9 

39.9? 

/7jrs 

.99 

/3.S-& 

3-9/. 

f*b" 

S3£/3 

9r,<f 

3s-  3 

/£-(/ 

<?h 

/S7  73- 

330 

/sr->~ 

S3H/S 

f%3 

3?,^ 

/Z9# 

-9/ 

7/,0:/ 

37/ 

/8c 

o^/f 

pes 

,*£* 

X/.93 

*&. 

20.9  7 

<Vs~7 

ss 

^3^/3 

fr.3 

33.  -5_ 

/3>3lC 

>$? 

77.7.9 

ZZZZ 

■9x 

66- 

536ft 

9f.c/ 

3  s-  3 

/£./ 

<9$ 

//■/7L 

PV7 

c7w 

76~" 

538  & 

/ob,  9 

3>^ 

70,  71 

-59 

9-77 

3S* 

*6S- 

77 

SVc  3 

/C7./ 

43 

9-93 

A  00 

9.93 

7S.3 

■  s~c 

SO 

S  9*3/3 

//*,-$ 

V^.^L 

Z-&7 

7.Q/ 

7.  ?5- 

7  7/ 

o^S     jy 

8/ 

S*r* 

//s~/ 

s~/.a 

$.37 

70  X 

93s- 

o&0 

'W 

^rS 

/5W/P 
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APPENDIX  B 
EDGE  EFFECTS 

In  the  set  of  experiments  involving  the  cylinders  con- 
structed from  Temsheet,  the  top  and  bottom  edges  of  the  test 
section  were  not  guarded  against  heat  losses.   An  analysis 
was  performed  and  it  was  found  that  edge  effects  did  not  in- 
fluence the  temperature  field  beyond  a  distance  of  approxi- 
mately 1/4  inch.   The  details  of  the  analysis  are  presented 
in  this  appendix. 


Top     — > 

Electrode 

y  \ 

\ 

L 

\ 

Bottom 

r.  lec  i-iroae-* 

w 

f 

The  axial  variation  in  the  temperature  field  due  to  end 
losses  was  estimated  by  treating  the  Temsheet  as  a  thin  fin. 
Symmetry  about  the  centerline  was  assumed.   Temperature 
variations  in  the  circumferential  direction  were  neglected. 
A  "worst  case"  situation  was  studied  by  assuming  the  temper- 
ature at  y=L  equaled  the  free  stream  air  temperature,  T^. 
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Due  to  symmetry,  the  center line  represented  an  adiabatic  sur- 
face.  The  governing  equation  for  this  situation  is: 

&\    +  J'  "  -  £-  (T-T  )  =  0     (1) 
dyz   k      kt     °° 

where 

T  =  Temperature  at  location  y 

q    =  Volumetric  heat  generation  rate  =  V2/2RLWt 

V  =  Voltage  impressed  across  the  electrodes 

R  =  Electrical  resistance  of  material  between  electrodes 

h  =  h  +h„  =  Surface  heat  transfer  coefficient 
c   r 

k  =  Thermal  conductivity  of  test  section  material 
t  =  Thickness  of  test  section  material 
Too  =  Air  stream  temperature 

The  boundary  conditions  on  the  problem  are: 
at  y  =  0,  dT/dy  =0  (2) 

at  y  =  L,  T  =  T„  (3) 

It  should  also  be  noted  that  the  surface  of  the  Temsheet 

that  was  in  contact  with  the  glass  wool  was  assumed  to  be 

perfectly  insulated.   The  solution  to  equation  (1)  that 

satisfies  boundary  conditions  (2)  and  (3)  is 

T-Too  =    {,         Cosh  (my)  \ 

T^Tco  Cosh(mL)  (4) 

where 

•  / 1 1 

(T-Too)    =   Centerlme   temperature   excess   =   q t 

h 

=1£ 


m  = 


Equation  (4)  was  used  to  determine  the  minimum  distance 
from  the  electrodes  that  data  could  be  taken  and  not  be 
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affected  to  any  substantial  degree  by  end  losses.   As  a  guide- 
line, equation  (4)  was  used  to  determine  the  y  location  for 
which  T-T  /T0-Too  =  0.99.   Using  representative  values  of: 
2L  =  6  inches 
h  =  5  BTU/hr-ft2-°F 
k  =  0.06  BTU/hr-ft-°F 
t  =  0.00325  ft 
This  distance  is  found  to  be: 
y  =  2.78  inches. 
It  was  concluded  that  if  the  electrodes  were  spaced  6 
inches  apart,  using  only  the  center  4  inches  would  yield 
temperature  information  that  was  essentially  uninfluenced 
by  edge  effects.   Figure  14  in  the  body  of  the  Thesis  sup- 
ports this  claim. 
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APPENDIX  C 
COMPARISON  OF  THERMOCOUPLES  AND  LIQUID  CRYSTALS 


In  the  set  of  experiments  that  employed  the  Nichrome 
wrapped  cylinder,  a  comparison  was  made  between  the  tempera- 
ture determined  with  the  thermocouples  and  the  temperature 
as  determined  visually  using  the  liquid  crystal  tapes.   This 
comparison  was  made  to  serve  as  a  check  on  the  accuracy  with 
which  one  can  determine  temperature  using  liquid  crystals. 
In  comparing  the  two  temperature  sensing  techniques  it  must 
be  noted  that  the  thermocouples  were  located  on  the  inside 
of  the  Nichrome  ribbon  whereas  the  liquid  crystals  were  lo- 
cated on  the  outside  of  the  tape.   As  such,  a  thermal  resis- 
tance exists  between  the  two  locations.   This  will  result 
in  a  temperature  drop.   An  elementary  heat  transfer  analysis 
on  the  composite  system  as  shown  sketched  below  yields  the 
following  expression  for  this  temperature  drop 


TT 


i 


Tape  Composite 


Nichrome i 

//////////tin  1 1  urn  i  mini) 


T 


k2LL  (1  L  >     +2ki(L'  J 
The  subscripts  "1"  refers  to  the  Nichrome  and  the  sub- 
script "2"  refers  to  the  tape.   The  volumetric  heat 
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generation  term,  q''',  is  related  to  the  voltage  drop  across 

the  ribbon  as  q   "   =  V2/RAt 

where 

R  =  Electrical  resistance  of  the  ribbon 

A  =  Total  surface  area  of  the  ribbon 

t  =  Thickness  of  the  ribbon 

V  =  Voltage 
The  following  were  taken  as  typical  thicknesses  and  properties 

t1    =  .00025  ft 

t2  =  .00075  ft 

kj  =  5.5  BTU/hr-ft-°F 

k2  ■  .12  BTU/hr-ft-°F 
Using  these  values  the  expression  for  AT  becomes 

AT  =  2.07  x  10_6q' ' '   [°F] 
where  q    ' '    is  expressed  in  BTU/hr-ft3 

The  above  expression  was  then  used  to  evaluate  the  re- 
sults of  the  test.   Several  items  where  noted.   The  thermo- 
couple readings  should  have  been  higher  than  those  of  the 
liquid  crystals.   This  was  not  always  the  case.   Careful 
analysis  of  the  surface  showed  that  the  thermocouples  were 
only  under  tape  R-45.   The  liquid  crystal  temperature  dis- 
tribution was  taken  from  several  tapes.   This  meant  that  the 
temperatures  compared  were  at  the  same  angular  location  but 
not  necessarily  the  same  vertical  location  on  the  cylinder 
surface. 

It  should  be  noted  that  the  location  of  the  thermocouples 
under  tape  R-45  was  easily  found.   The  thermocouples  acted 
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as  a  heat  sink  and  could  easily  be  found  during  warm-up  of 
the  cylinder.  One  showed  up  as  a  pin  point  size  cool  spot 
at  an  angular  location  of  90°. 

Examination  of  the  surface  of  the  cylinder  showed  that 
the  thermocouple  was  located  on  a  smooth  surface  whereas  the 
surface  under  the  tapes  directly  above  or  below  it  could 
have  irregular  shape  because  of  the  Nichrome  ribbon  irregu- 
larities.  The  surface  irregularities  produced  hot  or  cold 
spots  and  accounted  for  the  departure  from  expected  results. 

Since  thermocouple  position  on  tape  R-45  was  known  and 
the  surface  there  was  smooth,  it  was  decided  to  compare  tem- 
perature readings  at  that  point.   The  readings  compared 
within  the  calculated  differences.   Table  II  shows  this 
comparison. 

Table  II 

Comparison  of  Thermocouple  and 
Liquid  Crystal  Temperature  Readings 


Temperatures  (°F)         Temperature  Differences  (°F) 
Liquid  Crystals   Thermocouple   Maximum  Predicted  Actual 


108.9    (Red) 

109.8 

109.9    (Green) 

112.5 

112.1    (Blue) 

113.5 

2.7 

.9 

2.9 

2.6 

3.2 

1.4 
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APPENDIX  D 
UNCERTAINTY  ANALYSIS 

The  degree  of  uncertainty  for  the  final  results  was  cal- 
culated using  the  method  of  Kline  and  McClintock  described 
in  Ref.  20. 

The  measured  variables  which  were  the  origins  of  uncer- 
tainty were  voltage,  resistance,  surface  temperature,  angular 
location,  micromanometer  pressure  drop,  cylinder  physical 
measurements,  and  the  properties  of  air. 

The  error  in  determining  the  properties  of  air  was 
negligible. 

THE  HEAT  TRANSFER  COEFFICIENT 

The  heat  transfer  coefficient  was  calculated  from 

h  =  V2/R 

AAT 


(2aV)   +  fe2  +  (^AT)   +  (^A) 

h    II    V        R  AT        A 

For  the  following  experimental  run  at  a  Reynolds  number  of 

47,818 

V   =  27.6  ±  .3  VOLTS  (20  to  1) 

R  =  12.0  ±  .3  OHMS  (20  to  1) 

AT  =  39.7  ±  ,9°F  (20  to  1) 

A   =  .52  ±  .01  ft  (20  to  1) 
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6>  v12;  ^39.7'  \52; 


Substituting 

s>  ."VIZ 

h     \|v27. 

wh  =  .045 
~h 

This  represents  the  uncertainty  at  the  forward  stagnation 
without  regard  to  the  uncertainty  in  determining  the  angular 
location.   It  is  estimated  that  the  uncertainty  in  reading 
angular  location  in  this  region  would  be  ±5°  (20  to  1) . 

The  uncertainty  in  the  heat  transfer  coefficient,  at 
separation  was  approximately  4.5  percent  with  an  angular 
uncertainty  of  ±1°  (20  to  1) . 

The  uncertainty  in  the  wake  was  on  the  order  of  8  percent 
with  an  angular  uncertainty  of  ±5°  (20  to  1) . 

THE  NUSSELT  NUMBER 

The  Nusselt  number  was  calculated  from  the  equation 

Nu  =  hD 
K 

The  uncertainty  equation  would  be 


'  Wv,   2   ■    Un  2   ,    U)v  2 

( _£)   +  (  D)   +  (  K} 


^Nu  = 
Nu     y  '  h'      '  D'      v  K' 

Since  the  uncertainty  in  D  and  K  was  considered  negligible 


u 


Nu  =  wh 


Nu    h 
or  the  same  as  for  heat  transfer  coefficients. 
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THE  REYNOLDS  NUMBER 

Reynolds  number  was  calculated  from  the  equation 

Re  =  VD 
v 

and  its  uncertainty  is  the  same  as  free  stream  velocity,  V. 


-1 


2Ap 
P 


CO 


Re  = 


Re 


Vl 


(.5%>}2  =  .5("AP) 


Ap' 


Ap' 


Since  Ap  is  directly  related  to  manometer  reading  Ah 

"Re  =  .5  "Ah 
Re        h 

For  a  Reynolds  number  of  37,885 

mRe  =  .5  ..01   =  .025 
Re      \19J 

For  a  Reynolds  number  of  148,000 


0) 


Re  =  .5  ,.01.  = 
Re      l2.8J 


.002 


THE  FROESSLING  NUMBER 


The  Froessling  Number  was  calculated  from  Fr  =  — 


Then, 


to 


Fr  = 


Fr 


MNu 2  +   "Re  2 
1  Nu'     l2Re' 
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Substituting 

^Fr  =   l|(.045)2  +  (.017)2  =  .05 
Fr     ' 

This  is  for  the  Reynolds  number  47,818  used  previously.   It 

would  decrease  with  increasing  Reynolds  numbers. 
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